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Abstract. 

ON We study the formation and the evolution of small groups of galaxies 

1—1 using new Monte Carlo simulations. These are directly based on the 

+^ random walk approach to the statistics of condensations collapsing by 

gravitational instability, and the results are consistent with the Press & 
Schechter formula. We stress how observational features of groups, such 
as the galaxy membership, sensitively depend on the global dynamics 
of the universe. This is because the amount of aggregations between 
infalling galaxies is governed by the ratio of the group age to the time scale 
for dynamical friction; this ratio in turn strongly depends on cosmology. 
Thus in low-density universes - even if flat - the aggregations would lead 
often to a strongly dominant merger, contrary to the observations of 
t-h compact groups. 

o\ 

^ H 1. Introduction 

jH We submit that the galaxy population of small groups of galaxies (M ~ 10 13 M ) 

^ constitutes a sensitive probe of the global dynamics of the universe. 

To construct a realistic model for the evolution of such a population, we 
start from the results of Governato, Tozzi & Cavaliere (1996, hereafter GTC). 
There we showed with N-body simulations that the structure of small groups 
results from the interplay of three basic processes: collapse of the initial density 
perturbation, aggregations between the member galaxies, and secondary infall of 
additional galaxies. Our main result was that in a critical universe the "merging 
runaway" (see Cavaliere, Colafrancesco & Menci 1991) is over by the time when 
secondary infall is still on, so that the galaxy population of the groups is effec- 
tively rejuvenated. This appears to be a key to understand the very presence 
today of compact groups (Hickson 1982). 

On the other hand, in open universes the suppression of secondary infall 
leaves the dominant role to aggregations between galaxies, which lead to a single 
merger remnant (see also Tozzi, Governato & Cavaliere 1996). 
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Figure 1. a) Random walk in the S — S plane; b) the corresponding 
mass history in the M — t plane (arbitrary units). 

Here we intend to examine further the role of aggregations in small groups 
of galaxies for a wider range of cosmologies and of power spectra for the initial 
perturbations. For the sake of clarity we refer by "aggregation" to the coales- 
cence of the baryonic components of the galaxies, as distinct from "merging" of 
dark halos. 

2. Mass histories from random walking density perturbations 

We use an efficient way to compute a large number of mass histories of groups, 
constituted by Monte Carlo simulations. There is no unique way to generate 
such mass histories. Lacey & Cole (1993) proposed a Monte Carlo which "tends 
to overestimate halo ages" relative to N-body simulations. Kauffman & White 
(1993) constructed merger trees by extracting a posteriori the progenitor mass 
distribution from the analytical formula (Bond et al. 1991, Bower 1991). 

We propose a new version directly based on the random walk description 
of the statistics of the perturbation field, as proposed by Bond et al. (1991). 
We start directly from the random walks (trajectories) in the plane S — S, S 
being the density contrast at a resolution S(M) given by the variance of the 
linear field in the form S(M) = a 2 (M). The increment in S after a step dS is 
a stochastic variable extracted from a Gaussian with variance dS. Once S first 
exceeds a given threshold, an object appears at the "time" S(z) at the variance 
S(M) (see fig la,b). For every merging event, i.e. for every discrete variation 
in the mass of the group, a new branch appears in order to conserve mass, so 
that a complete merger tree can be constructed. Our merger trees are tested for 
consistency with the analytic mass distribution of the progenitors (Bond et al. 
1991, Bower 1991). 

We stress that, at variance with previous simulations, our Monte Carlo is 
based on a true Gaussian process. Moreover, merging events univocally corre- 
spond to a discontinuity in the mass history M(z) (i.e., they do not depend on 
dS, see fig la). We also stress that this procedure is consistent with the Press & 
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Figure 2. Average mass histories of the main progenitor. Solid line: 
TCDM; dotted line: OCDM; dashed line: ACDM. 

Schechter (1974) formula; in particular, the same dispersion of M — z is obtained 
with the only requirement of mass conservation in the merger tree. 

Next we present an application to the evolution of groups. For simplicity, 
we follow only the most massive branch of each merger tree, i.e. the "main 
progenitor" . 

3. Histories of small groups of galaxies 

For the power spectrum of the initial perturbations, we select three represen- 
tative CDM models: the tilted CDM model in the critical universe (TCDM), 
the open CDM with Q = 0.2 (OCDM), and finally the flat-low density CDM 
model (ACDM) with Q = 0.2 and A = 0.8 (Bardeen et al. 1986, White et al. 
1996). The normalization of the power spectra, expressed in terms of a s (the 
rms value of the linear fluctuations on scales 8/i _1 Mpc) are calibrated using 
the local observations. Among them most solid are the local abundance and the 
temperature distribution of clusters (White, Efstathiou & Frenk 1993, Henry 
& Arnaud 1992). These values are: a 8 = 0.6 for TCDM, a 8 = 1 for OCDM 
and <7 8 = 1.4 for ACDM. These are in agreement with the normalization given 
by COBE for such spectra, within the statistical uncertainties (Bunn & White 
1995). Moreover, we take H = 50 km sec _1 Mpc" 1 for TCDM and H = 65 
km sec - 1 Mpc" 1 for OCDM and ACDM. 

Following the most massive branch of each tree, we compute the average 
history of the main progenitor (see fig 2). Next, we compute the distribution 
of formation epochs Zf, defined as the redshift at which the main progenitor 
reaches one half of the final mass; equivalently, we present the distribution of 
the ages of the groups counted from Zj (see fig 3). 

In the critical cosmology present-day halos were assembled only recently, 
gaining most of their mass by infall of dark matter halos in the fairly recent 
past. The formation epochs are peaked around Zj ~ 0.5, and about the 85 % 
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Figure 3. a) Distribution of the formation redshifts Zj. b) Ages of 
the main progenitor in Gyr. 

of the present-day population is younger than z ~ 1, i.e., younger than about 8 
Gyr. 

In OCDM the majority of the merging events, and so most of the mass 
growth of the halos, occur at high redshifts. In this case the values of Zj are 
widely distributed around Zj ~ 1.5; more than 80% of the groups are older than 
8 Gyr. In ACDM the distribution is intermediate, with average Zj ~ 1. 

It is then clear that in the critical case the galaxy population of the groups 
is rejuvenated at recent epochs by the galaxy included in the infalling matter. 
There is little time for these galaxies to suffer extensive dynamical friction, and 
so low probability to aggregate with the merger remnant. 

In OCDM and in ACDM, on the other hand, the galaxy population is 
formed at earlier epochs. So the majority of the galaxies spend most of their life 
in the group with enhanced probability to suffer aggregations. 

In fact, the time scale which governs the onset of aggregations inside groups 
is the dynamical friction time scale Tj. The classic formula for Tj (see Binney 
& Tremaine 1987), which is in agreement with the behavior seen in our N- 
body simulations (GTC), is valid when the infalling galaxies have masses much 
smaller than the receiving group, a condition satisfied by the vast majority of 
the infalling halos. If we assume (following Lacey & Cole 1993) that on infalling 
each halo is stripped of its dark matter component, we obtain a lower bound to 
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Figure 4. Fraction of present-day groups with n surviving galaxies 
in different cosmological models (the lines are the same as in fig 2): a) 
for elongated orbits of the infallen galaxies; b) for circular orbits. 



the number of accreted baryonic cores, since the infalling halos can host more 
than one baryonic core (with less dynamical friction). 

To give an example, we consider halos infalling into the main progenitor for 
z < Zf, and we calculate Tj for halos with M > 10 12 M Q . One such event is 
interpreted as one bright galaxy infalling into the group, assuming that each halo 
contains a galaxy with M g = /jM, where f b = Q b /Q is the baryonic fraction 
(here we assume Q b = 0.013 h~ 2 consistent with nucleosynthesis, see Walker et 
al. 1991). The ratio of the group age to Tj sets the amount of aggregations of 
galaxies with the central merger remnant. In fact, if Tj is smaller than the look- 
back time to the infall event, the galaxy may be assumed to suffer aggregation 
with the central merger remnant. In this way we obtain a lower bound to the 
number of survivors. The result depends on the elongation of the orbits (see 
GTC). For elongated orbits and TCDM, we see that more than 75% of groups 
contains at least a bright surviving galaxy, while in OCDM less than 30% of 
groups contain a bright surviving galaxy. In ACDM the fraction is raised to 
37%. For circular orbits, which minimize the dynamical friction, there is still 
difference between the models (see fig 4). In this sense the groups behaves like 
cosmological clocks (see fig 5). 
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Figure 5. Small groups as £7 — clocks. The arrows indicate cosmic 
time measured in units of Tf (represented with the dashed area), to 
visualize that with TCDM, only a few units of Tf have elapsed since 
z = Zf] with OCDM instead, the present-day groups are many units 
of Tf old, and the aggregation process goes to its completion. 

4. Conclusions 

We have studied the link between the structure of small groups of galaxies and 
cosmology. We stress that such link stems from the interplay of two basic time 
scales: the age of the small groups (set by the cosmological model) and Tf for the 
infalling galaxies. These time scales are comparable for small groups of galaxies 
with M ~ 10 13 M Q . 

The groups constitute sensitive cosmological clocks because their structure 
is affected twice by the cosmological model: as for the distribution of formation 
epochs, and as for the degree of aggregation between infallen galaxies. 

To examine these processes, we have presented a new, efficient Monte Carlo 
scheme to extract mass history of halos directly from the statistics underlying 
the P&S formula. 

Our preliminary results predict that in the critical universe infallen galaxies 
have a good probability to escape aggregations down to the present time, and 
so constitute true groups of few to several members. In low-density universes, 
in which the majority of small groups form earlier and afterwards experience 
little secondary infall, the galaxy population is not rejuvenated; aggregations 
constitute the important dynamical process which goes to its completion and 
forms a single remnant or a dominant merger surrounded by satellites. Flat low- 
density universes constitute intermediate cases, but closer to the open version. 

We thank N. Menci for many useful discussions, and P. Padovani for his 
help. We acknowledge research grants from MURST and ASI. 
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